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EMISSION OF SOUND FROM TURBULENCE CONVECTED BY A PARALLEL 
FLOW IN THE PRESENCE OF SOLID BOUNDARIES 
by Marvin E. Goldstein and Burt M. Rosenbaum 
Lewis Research Center 

SUMMARY 

A theoretical description is given of the sound emitted from an arbitrary point in a 
parallel or nearly parallel turbulent shear flow confined to a region near solid bound- 
aries. The analysis begins with Lighthill's formulation of aerodynamic noise and as- 
sumes that the turbulence is axisymmetric. It leads to equations which predict the 
directivity patterns of the spectral density of the intensity of the sound emitted from a 
unit volume of turbulence located at an arbitrary point within the flow in terms of certain 
correlation- length scales of the turbulence. In order to apply these equations to predict 
the sound field emitted from a complete flow, it would be necessary to integrate the 
results over the flow field. Specific results are obtained for the sound emitted from an 
arbitrary point in a turbulent flow within a semi- infinite open-ended duct. An experi- 
ment which could be used to verify these results is proposed. 


INTRODUCTION 

Since STOL transport aircraft are intended to operate in densely populated areas, 
it is necessary to give considerable emphasis to the problem of noise reduction in the 
development of these aircraft. Most STOL aircraft depend for their performance on 
lift augmentation devices which involve high-velocity turbulent flows along solid sur- 
faces. It has been found (ref. 1) that the presence of such surfaces can result in a con- 
siderable increase in noise production. 

In this report we shall attempt to obtain some insight into the effect of solid bound- 
aries on the process of aerodynamic sound generation and at the same time make some 
progress toward developing techniques for predicting such effects analytically. To this 
end we shall develop theoretical formulas for predicting the sound emitted from a unit 
volume of turbulence located at an arbitrary point in a nearly parallel turbulent shear 



flow confined to a region near one or more solid boundaries. Typical configurations to 
which the analysis applies are shown in figures 1(a) and (b). It should be emphasized 
that the analysis, by itself, only relates the sound field generated by the turbulence at 
an arbitrary point in the flow to the turbulence correlation length at that point. In order 
to use this analysis for actually predicting the sound field from a given flow, it is neces- 
sary to integrate the formulas over the entire region of turbulence generating this sound. 
Since the turbulence correlation length appears under the integral sign in this integra- 
tion, it is necessary to know this correlation length at every point of the flow in order to 
predict the sound field. However, determining the distribution of this correlation length 
is in itself a difficult problem, and it will usually be necessary to rely heavily upon ex- 
periment for this purpose. Nevertheless, the formulas obtained can be used (even when 
the distribution of the correlation length is unknown) for studying the effect of the solid 
boundaries on the sound field from the turbulence located at any individual point within 
the flow. 

The analysis begins with Lighthill’s model for the sound emission from turbulent 
shear flows (ref. 2). The turbulence is assumed to be locally homogeneous and axisym- 
metric; that is, the statistical properties of the turbulence are independent of rotations 
about the direction of mean flow. However, they can be different in the direction of 
mean flow and the directions perpendicular to the mean flow. The analysis is carried 
through by introducing an appropriate Green's function for an arbitrary geometry and 
then using the symmetry properties of the turbulence to transform the result into a 
manageable form. It is then shown how the reciprocity principle can be used to simplify 
the problem of finding the Green's function for the case where the observer is in the far 
field. The final result of this analysis is a set of formulas for the spectral density of 
the far- field intensity of the sound emitted from a unit volume of turbulence located at 
an arbitrary point in the flow. These formulas can be interpreted as the sound field 
from a particular sum of statistically independent convected quadruples (with certain 
coupling terms) whose strengths are determined by the statistical properties of the tur- 
bulence . 

A particular lift augmentation device which is being considered for use in STOL 
aircraft is the augmentor wing ejector flap shown in figure 2. This flap is fed (blown) 
through the duct in the wing, as shown in the figure. The high-speed jet issuing from 
the wing into the flap (at point 1 in fig. 2) is a potentially important source of noise with- 
in the flap. The effect of the flap boundaries on the sound emission from this noise 
source can be important. 

Motivated by this problem, we applied the general analysis previously described to 
predict the noise emitted from a unit volume of turbulence located at an arbitrary point 
within an open-ended semi- infinite duct. The Green's function needed for the solution 
to this problem is obtained by transforming the boundary- value problem which deter- 
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mines this function into one whose solution is known. The turbulence properties which 
are used in this example are taken from observed results for the mixing region of jets. 
This is done in order to relate the analysis to the noise emitted from the turbulence in 
the high-speed jet from the duct in the wing into the flap. It must be emphasized, how- 
ever, that in order to use this analysis to predict the noise from the jet mixing region 
in a complete augmentor wing, it is necessary to know the distribution of correlation- 
length scales within this region and to integrate the results over the entire region of 
turbulence. In addition, in any real augmentor wing there will be additional noise 
sources such as the mixing region between the flow leaving the flap and the exterior flow 
about the wing and the possible noise source due to the flow entering the flap from above 
the wing. However, the general analysis is also capable of dealing with many of these 
additional noise sources, provided enough details of the turbulent flow are known. Thus, 
the present analysis, by itself, certainly does not give a complete solution to the prob- 
lem of predicting the noise from an augmentor wing jet flap. Nevertheless, it does pro- 
vide a necessary first step toward this goal. In addition, numerical results are pre- 
sented for the directivity patterns of the spectral density of the far- field intensity 
emitted from a unit volume of turbulence located at various points within the duct. 

These results show how the sound emission processes and the geometric effects of the 
duct interact. It also shows how the elementary sources of sound, that is, the individual 
turbulent eddies, contribute to the sound field. An experiment which should permit 
comparison with the computed results is discussed. 

In order to help differentiate which aspects of the sound field are caused by the con- 
vected quadrupole nature of the turbulent sound sources and which aspects are geometric 
effects of the solid boundaries, results are also presented for the sound emission from 
a simple monopole source located within the duct. 


FORMULATION 

We shall consider the sound emission from a nearly parallel turbulent shear flow 
confined to a region near solid boundaries. Typical configurations are shown in figure 1. 
The mean velocity is taken to be in the yj- direction, and it will be assumed that it var- 
ies only in the y 2 -direction. We begin with Lighthill's equation (ref. 2) for aerodynamic 
noise generation 


V 2 p - 


d 2 p _ _ 


2 ,.2 

c o 9t 


.2 3y.3y. - 


(p Vj + p ij - vy 


a) 
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where 



a 2 


(All symbols are defined in appendix E. ) Since aerodynamic noise generation is usually 
of interest only at very high Reynolds numbers, Lighthill (ref. 2) neglected viscous ef- 
fects. Following Ffowcs Williams (ref. 3) we shall suppose that viscous effects can 
also be neglected even when solid boundaries are present. We can therefore set p^ 
equal to pS^. We shall also neglect (ref. 2) the effects of heat conduction so that 

P - P 0 = c 2 (p - p Q ) 


Finally, we suppose that (ref. 2) 


P v i v j 


0 V - V * 
1 ] 


This assumption amounts to neglecting the direct effects of convection and refraction of 
the emitted sound by the mean flow and of the scattering of the emitted sound by the tur- 
bulence. 

When these approximations are introduced into equation (1) we obtain 


V 2 p - 


1 


0 2 p 


„2 -.2 
c o St 


3 % 

dy^Yj 


( 2 ) 


where 


T ij * Wj <3) 


In order to ensure convergence of the integrals to be encountered, we shall suppose that 
t. • = 0 for t > T, where T is some large time interval which will be put equal to in- 
finity at the end of the analysis. 

We shall suppose that the solid boundaries are rigid so that the normal velocity 

v = v-n- (where n- is the i^ 1 component of the unit normal directed outward from the 
n i i v 1 

turbulent region at the solid boundaries) vanishes at the boundaries. With viscous ef- 
fects neglected, the equation for the normal component of momentum at the solid bound- 
aries becomes 
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V V u\. ~ 

p — - + pv-n. — - + n. — = 0 

at J 1 a Yj 1 3 Yi 

— 2 

The first term is zero. The middle term can be shown to be equal to ( +pv / R), where 
R is the radius of curvature of the boundary surface in the direction of the motion and 
where the negative sign holds when the boundary surface is concave inward along this 
direction and the positive sign holds when it is convex inward. Thus, we get 

n • Vp = ± 

R 


We shall assume that the surfaces adjacent to the regions of turbulent shear flow 
are relatively flat; that is, R >> h, where h is a characteristic length representative 
of the distances between boundary surfaces (see fig. 2). Then pV /R ^ 0 for surfaces 
adjacent to the flow. We shall also suppose that, at the surfaces outside the regions of 
turbulent shear flow, the acoustic approximation applies; that is, the terms involving 
the square of the velocity are negligible. Hence, 


n. -^2. = h • Vp « 0 


(on solid boundaries) 


(4) 


that is, the normal derivative of the pressure essentially vanishes at the solid bound- 
aries. 

It is convenient to work in terms of the Fourier transforms (with respect to time) 
of the variables. Hence, we put 


P = 


1 




e iwt (p - P Q ) dt 


(5) 




( 6 ) 


Then equation (2) becomes 


V 2 P + k 2 P = - 


i!Zii 

ay^yj 


( 7 ) 
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where 


k--H 


is the wave number. Then the boundary condition (4) becomes 

n • VP = 0 (on solid boundaries) (9 

The boundary condition at infinity is that P represents an outgoing wave. 

The solution to this problem (eq. (7) with boundary condition (9)) can be written in 
terms of the Green's function G, whose normal derivative vanishes on the solid bound- 
ary. Thus, 


/ 


/ a 2 T. • 

P(x ) = / G(x I y ) dy 




where V denotes the volume outside the solid boundaries 


V 2 + k 2 )G - -6(x - y ) 


n • V_^G = 0 (on solid boundaries) 
x 


where vi = 3 2 /0x^9x^ and G represents an outgoing wave at infinity. Differentiating 


by parts shows that 


(GT.J - j- (tj. -i- gV -L (V- j- g\ = G-^ii - T tj 

3 3Yi \ 3 3yj f y 1] a yi j dy^y. ] 0y i 3y j 

Hence, the divergence theorem shows that equation (10) can be written as 

P(x) = f T.-J-SLdy + f n.(-2- GT..^ dg- f n-T^-^-d^- f n-T^ dg 

Jy Jg Vi 7 Jg J 9y i Jg J ]9y i 
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where £ denotes the solid boundaries. But since 


t « =-p r ei “ tp o v i v j m 

V 00 


the requirement that the normal velocity vanish at the solid boundaries shows that the 
last two surface integrals on the right side of equation (12) vanish. And since the cur- 
vature of the solid surfaces adjacent to the region of turbulence (i. e., where T.j is 
nonnegligible) vanishes in the flow direction, the rate of change in the normal vector in 
the direction of flow (vj/| v | )(3/Syj)n is perpendicular to the flow direction v/|v | . 
Hence, it follows that 



(-i-L 

\ l v I 9 yj 



= o 


so that 


Vj 



and 


v,— (Vjjij) = n ( v. — 

’ 3y j J a yj 


Therefore, 


n,— GT„ =— !- / e‘ 

3y i 1 Vi ijL 

el 


n-Vj— (Gv.) + v.G— (n.v.) 

ay j 3 3 ^j _ 


dt 


v n i^ + Gv,-Lv n 
_ 3y i ’ 0y i J 


dt 


But since the second derivative is along the surface and v n is identically zero on the 
surface, the remaining surface integral on the right side of equation (12) also vanishes. 
Hence, equation (12) becomes 
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(13) 


P(x ) = 



3 2 G 

Sy^yj 


dy 


Although the integral is formally over all space, it need only be carried out over the 
source region. 

The spectral density of the intensity in the far field is 


I 



= P(x)| 2 

!T v» 


where T is some large period of time such that t.j can be considered to be zero out- 
side the time interval (-T, T), which is of duration 2T. At the end of the analysis, the 
limit will be taken as T — 00 . 


I (x) =■ 

or ’ 


1_ f f T ii (f,)T k 

o p o Jy Jy 2T 


y”) 


5 2 G(x y T ) 9 2 G* 


(x 


3y-9yj 3y!'8y!’ 


2—1 dy' dy" 


(14) 


Now upon using the fact that the Fourier transform of a convolution is proportional to 
the product of the Fourier transforms of its components and the fact that the Fourier 
transform of f(-t) is the complex conjugate of the Fourier transform of f(t), provided 
f(t) is real, it follows from equation (6) that 


T i]ty , r r kz<y > '> 

2T 



e" iwT v j v j v k v z’ 


(15) 


where 


v i v j v k v Z 



v^y’jtjv^y'.tjv^y",! + r)v z (y",t + r) dt 


is the fourth- order two- point two-time velocity correlation tensor. 


(16) 
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TURBULENCE CORRELATION FUNCTIONS 


The following development closely parallels that given in reference 4. Since the 
mean velocity U(y 2 ) is assumed to be in the y^-direction, the velocity v^ can be ex- 
pressed in the form 

v^y , t) = 6 il U(y 2 ) + Ujty ,t) 

where u. is the turbulent (fluctuating) component of the velocity. 

Upon introducing this into equation (16) we find that 

v i v j ^(yS^r) + i (17) 


where we have put 


rf = y " - y ' = 


(18) 


* o’VijW + + o ,2 o"Vi,« u »n +u i u j u k u F 


r) ■ u i ,u j u i' u i' - ujuj ujj'u;’ + 2U'(6 li u j ’u i -u i ' + 


- U’U"(6 li s lk ujuj , + 5^^ + (is) 

/J «L[^ ,,7<1) ’ T ] s u i u j u k u i’ - vj “k“? * 2 u "( 6 ik u i u j u i' + htm 


+ u ' u "(6iiSi k u:u;' + + 2S U « al «J»i’) (20) 


The double primes indicate that the quantities are to be evaluated at y" = y' + if and 
t + r, while the primed quantities are evaluated at y" and t. Notice that (y',rf) 

is independent of r. Hence, 


~ h e "'“ T *Si W <?> ’ 7) dT = 5M(! mi 


where 6(w) is the delta function. 
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At zero frequency the inhomogeneous Helmholtz equation (11a) reduces to Poisson’s 
equation. For |x| >> |y| (i.e., in the far field), the delta function vanishes, and the 
Green's function of the resulting Laplace's equation is independent of the source location 
y. Hence, 


— = 0 at u> = 0 
3Yi 


But this implies that 


6(o>)— - 0 

ayi 

It therefore follows that, when equation (17) is substituted into equation (15) and the re- 
sult is substituted into equation (14), the term ^h(P'>v) makes no contribution to the 
integral. Hence, upon making this substitution and changing the variables of integration 
from p and y" to y’ and if for the integration over and from p and y” 

to y" and if^ for the integration over we find that 



where V 2 denotes the appropriate volume for the double integration over the variables 
y’ and if (which is the same as for the integration over the variables y", if^). But 
since the variables y" and rj^ in the second integral are variables of integration, we 
can replace them by y’ and if, respectively; and since i, j, k, and l are dummy 
indices, we can interchange the indices i and k and the indices j and l without 
changing the value of this integral. Thus, the second integral in equation (21) becomes, 
after changing the variable t to - t, 


Pq 

4?rc o 



9 2 G* (x j p ) 9 2 G(x|f + ip/g(2) _ T)e icoT 

sy:a y : ay^ay^ 


dr dif dy ’ 


( 22 ) 
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Now, if it is assumed that the turbulence is stationary in time, the correlations must be 
invariant under a translation in time. Hence, we can replace t in \rf,-T) by 

t + r. It can now be seen that replacing y" by y\ rf^ by if, and t by -t inequa- 
tion (20) merely interchanges the primed and unprimed quantities. Hence, it follows 
from equations (19) and (20) that 

And upon inserting this into equation (22) we find (since the actual limits of integration 
on the variables y" and if^' are the same as those for y' and if^\ respectively) 
that the second term on the right side of equation (21) is the complex conjugate of the 
first. Hence, 



EXPANSION FOR LOW EDDY MACH NUMBER - CONVECTIVE AMPLIFICATION 

We first introduce the new variable 


L y 2 + y 2 

V 1 ’ 2 ’ 2 y 2 



where i denotes the unit vector in the y^-direction. Then 


y* 



(ii?j - v) 


Hence, if in equation (23) we put 

i&> *> t)= ^ijki ( y > ^ + T > 


where 
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and then change the variables of integration from y’ and if to y and if in perform- 
ing the integration over /^j^, we obtain 



where the integral involving was obtained simply by changing the name of the 

dummy variable of integration from y' to y. Thus, 


r _ r — — — — = — 1 (26) 

- 2U(y 2 )|5 1 iUj(y > t ) u k (y + if, t + T)Uj(y + rf, t + t) + 6^(7, t)u k (y + i 7 ,t+r)u z (y +if,t+T)J 

The symbol £/ 2 always denotes a given volume of space even though the specific limits 
for the transformed variables can be different. 

If the mean flow were zero, the correlation functions , t?, t) and 

would decay rapidly as functions of 77 relative to a characteristic length of the region of 
flow and the characteristic wavelength associated with the eddy decay time; and it would 
be possible (as will be seen) to neglect variations in G over the range of 77 where (l ^ 
and are nonzero. However, in the case where there is a mean flow (see ref. 5 ), 

the correlations at points along the direction of flow can extend over long distances com- 
pared with the wavelength associated with the eddy decay time since the eddies are moving 
with the flow and their apparent decay time is quite short. In order to allow for this, we 
introduce a set of coordinates which move with the turbulent eddies. Thus, we introduce 
the new variable 



where U c , which is in the yj direction, is the convection velocity of the eddy. The 
moving- axis correlation functions are now defined by 


R W^’ T) ^W^’ T) 

R iW^’ t) 

Then upon introducing the change of variable (27) into equation (25) we get 


( 28 ) 


y*)=— Re 

2 ire 


\fM 


2 o(x|y + iUi-ir) 3 2 C*(x|y + iUi + if + U c T) + 




ay k 3y. 


Rjjk;(y > ? . T )e dr dy d£ 


P f C£°l iii 

L Jyz) Sy i Sy j 


-.sVfxly + £~ + U r) _ _ _ _ 

^ R ijk£^ y > 4 > T > e ' dT dy d£ 


5y k 5y £ 


(29) 


Now let l denote a typical correlation length (eddy size) of the turbulence, let 
denote a typical decay time of the turbulence in the moving frame, and let be the 
Mach number characteristic of the turbulence. It is shown in reference 6 (see sec. 3, 
eq. (17)) that 


In turbulent shear flows the characteristic turbulent velocity is about l/10 th of the mean 
flow velocity. Hence, for subsonic flows 


T f c o 


« 1 


Thus, the eddy size is very small compared with the wavelength A f associated with the 
eddy decay time. For pipe flows and internal flows in general, the eddy size is very 
small compared with a typical dimension h of the boundaries (see fig. 1). Hence, we 
shall suppose that 


l« h 

The two length scales which appear in the Green's function G (when the observation 
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point is in the far field) are and h. Thus, it is expected that the Green's functions 

G^x |y Uf - fj, G* (x |y + f + U c r), and G*^x |y +-^ 1 ^ +-^ f + U c rj will change 

only slightly when £ varies over the range 0 ^ | £ | ^ l . But by definition 
R^jk^y, h , t) ~ 0 and R^j^y > K , r) « 0 for |i; | > l, so that we can neglect variations 
in the Green's function insofar as the integration with respect to i; is concerned. In 
addition, we can carry out the integrations with respect to £ over all space and the 
integrations with respect to y over the volume V exterior to the solid boundaries. 
When the approximations are introduced into equation (29), we obtain 



,-lCOT 


3 2 G*(x |y + ” 




—f [ R tVy-«' 


t) 


+ R ijkz ( y^’ T) ] d ^ dTd y {30) 

where the omission of a symbol for the limits of integration denotes that the integral is 
to be carried out over all space. 

We let I (x |y ) denote the spectral density of the intensity at the point x due to 
the sound emitted from a unit volume at y (ref. 7). Then 


d y 

and it follows from equation (27) that 

V* |y ) = 2^7 Gij( ^ ) / e " 1WTG ki ( * + Tl c T) /[ R ijkz(y >*> t) + R ijki ( y - t) ] d ^ dT 

(31) 


where 


) ■ 


s 2G „(?|y) 


Now it follows from equations (19) and (28) that 
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/ Kjk i - * . ■ r) ♦ r? j kZ (y.r, t)] dr = /s ljW <y , r , t) df 


/ u (y 2 - y) u (y2 + yj[ 6 ii 6 ik 6 J z(y>«‘- T ) + 2 Vik K i*<y-M df 


2U(y 2 )/[fl 11 B j( kZ (y , f . T) + «ijS ilti {y , F .■ T)] df 


(32) 


where 




Rij(y,r, t) s uju." 


Ri,j k (y,r, T ) s u^y^Ujfy + if,t + r)u k (y + ff,t + r) 


(33) 

(34) 

(35) 


The first two integrals on the right side of equation (32) represent the self noise and the 
shear noise, respectively (refs. 8 and 9). The last integral in equation (32) represents a 
coupling between the shear noise and the self noise. 


LOCALLY HOMOGENEOUS TURBULENCE - VANISHING OF COUPLING TERM 

It is shown in reference 4, under the relatively mild restriction that the turbulence 
is approximately locally homogeneous and incompressible, that the coupling term be- 
tween shear noise and self noise in equation (32) vanishes. Hence, after the names of 
the dummy indices are changed, equation (31) can be written as 


V x 


-I-. Pc 


y) = 




2n e 


Im> /; 


- 10 )T 




+ V) 


/Vz (f ’ r ’ T)d « dT ] 


+ 



re o 


G ij(y ) 



e’ lwT G 


11 


( y + u cT ) J v(y 2 +hy(y 2 .hy ji( yj t 


r)d t, dT 


(36) 
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The first term in this equation is the self- noise term and the second is the shear- noise 
term. 

Equation (36) cannot be further simplified without making additional assumptions 
about the turbulence. 

REPRESENTATION OF FOURTH-ORDER CORRELATIONS 
IN TERMS OF SECOND-ORDER CORRELATIONS 

It is argued by Batchelor (ref. 10, section 8.2) that the part of the joint probability 
of the velocity (with zero time delay) associated with the energy- bearing eddies is ap- 
proximately normal, at least insofar as the velocities at two points are concerned. This 
approximation is better for some purposes than for others. For example, this approxi- 
mation gives reasonably accurate predictions about the relation between the second- and 
fourth- order correlations (ref. 10, p. 176). This relation is found to be (ref. 10, 
eq. 8. 3. 11) 


i u j u k Z " i j k Z i k j Z + u i u Z u j u k 


at t = 0 


(37) 


But by extending the reasoning used by Batchelor in section 8.2, we can argue that, when 
the velocity correlations are separated in time as well as in space, their correlation will 
be subject to even more random influences from the neighboring flow than when they 
occur at the same time. These influences will, according to the central limit theorem, 
tend to make the joint probability distribution more normal. Hence, we expect equa- 
tion (37) to be even more nearly valid when t * 0. Then in view of equations (33) and 
(34) we can now write 


^ijkZ ^ ’ T ) - K ik^jZ + ^iZ^jk 


(38) 


Inserting this into equation (36) and changing the dummy indices gives 


yx|y) =— Re 
7rc 0 


* V>/ S ik«j l & dT 


2 P, 


7TC, 


U fie 


G^y) J e‘ iwT G* z (y + U c t) J U'U"R jZ df dr 


(39) 
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where U' s u 


^2 


and U" = U| 


^2 


AX I SYMMETRIC TURBULENCE -QUADRUPOLE CORRELATIONS 

It can be seen from the measurements presented in references 11 and 12 that, at 
least for jet mixing regions, it is reasonable to assume that the moving-axis correlation 
tensor R.j(y’, t) is a symmetric function of r. Thus, figure 27 of reference 11 shows 

that Rji(y t) peaks at r = 0, and figure 13 of reference 12 shows that this function 
is fairly symmetrical about r = 0. (Notice that in this figure the lines £ ^constant are 
parallel to the dashed line. ) Hence, we assume that R^ is an even function of r. We 
shall also assume that the turbulence is axisymmetric about the direction of mean flow. 
Then it is shown in reference 4 that R^ is an axisymmetric tensor of the form 

Rij(y > t, T) = M j + B6.. + C5 li 5 lj + + 6^.) (40) 


t 2 . t 2 . t 2 

^1 + h + ^3 

and are even functions of and where D is a function of y, r, and tj and an odd 
function of £ ^ • 

In order to evaluate the first term on the right side of equation (39), observe that 


where the coefficients A, B, and C are functions of y , r, and £ = 


V 


/Mu dr -/K*k + B5 ik + c 6 n 5 ik - D < 6 n«k - 6 ik«i>] 

X [ A {j{; + BCjj + CSjjS,, + D(6 lj{ ; + 6 u |j)]dT (41) 

All the products of the coefficients A, B, etc. which occur in the integrand are even 
functions of and £ 3 and either even or odd functions of £^. (In fact, they are all 
even except those terms which consist of products of D with either A, B, or C. ) But 
since any integral of an odd function of 4 . for i = 1, 2, 3 is equal to zero, it can be 
seen by inspection that only terms which have two pairs of equal indices can contribute 
to the integral. Hence, the integrals must be of the form 


I 


W d « = a < 


6 . . 6 , , 

lk i] k l 


+ b. . 6 ., 6. 7 
X] lk ]Z 


c. . 6 . , 6 ., 
i] iZ jk 


(no siim on i, j, k, or Z) (42) 
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Therefore, upon noting that = Gj. and interchanging names of dummy indices we 
find that 

Gjjty )Cj,(y + u c t) f R^Rjj df - ^ k G u (y )G^(y + u c t) + (b,, ♦ ^GyfylG*!? * u c r> 


It is convenient to define 

G i-(y )s G ij ( y + u c T ) (44) 

T 

This shows that the coefficient a ik of G^Gj^ is the sum of all terms in 

J * R ik R. ^ d|* with i = j and k - Z. Similarly, the coefficient b- + c^ of G^G? is 

the sum of all terms in J df with k and 1 equal to either i or j. There is 

a certain ambiguity in the definition of a ik and b^ + c- since those terms in 

J* R^Hj^ df* with all four indices equal can be included in either or these terms. But 

once a., is defined, b. + c. - is unique. These results are used in appendix A to show 
LK 13 1J 


= G ii G H/ s n d f 


( G 22 G 33 + G 33 G 22' _/"'^23 " ^12* d ^ + G 11 G 11 f ^ 11 ' ^12 ! d ^* 
2 ^ G 12 G 12 + °13 G ?3> /<«? 2 + K 11 S 22> ^ 


+ ,G 22 G 22 + G 33 G 33^ f ^22 " ^12* 

* 2G 23 G 23 /<* 22 - S 23> d «~ 

Applying the reciprocity principle to (11a) shows that for x * y 


G ii = " k G 
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Hence, the first term in equation (39) can be expressed in terms of the relation 

y^ik^j l d ^* = k4 ° GTs “ ( G 22 G 33 + G 33 G 22^23 “ ^22^ + G ij G ^j ( 45 ) 

where 




(46) 


and 


«11 f < K 11 - *?2> <f 


Q 12 " Q 13 


! / (B?2 + S ll«22> 


^22 “ ^33 ~ f & 22 " ® 12 ^ d ^ 


■/' 


q 23 = / (it*, - r|„) df 


22 " 23' 




(47) 


s-fn \ 2 dr 


(48) 


These five integrals depend on the four correlations R-q, **22’ **12’ and **23* 
Turning to the second integral in equation (39), we find from equation (40) that 


°ij<y> G l;<y + V> / u ’«" a i2 d «~ = G ij G u / u ' u "[ A «i«2 + B5 j; + C6 lj6l; 

+ D({ j 6 u + | ; 6 ij)] d r (49) 

It is shown in appendix B that this equation can be simplified to obtain 

G lj G L / u ’ u ” S j2 Of = G lj G lj / v ' v " S jj <«■ (50) 

It now follows from equations (39), (45), and (50) that 
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( 51 ) 


I w (*|y ) = I { „ ( *|y ) se + I w ( 5 r ly ) sh 

where the self- noise term is 

yx|y) ^f"k 4 C / e' iwT G T SdT 

7TC L 


-/ 


(G 


r T 

22 u 33 


G 33 G 22 


)(Q 23 - Q 3 2^ e 


■1WT 


dr 


G ij 


f 


-icoiy/IV. . 

CyQij dr 


and the shear- noise term is 


2 Pr 


V* ly) S h=— 

TIC ^ 




U'U’'R.. df dr 


(52) 


(53) 


The last term in the self- noise term (52) is a sum of the intensities from independent 
quadrupoles. The first term is the emission from a monopole source (weighted accord- 
ing to wave number), and the second term is a coupling term between the two longitudi- 
nal quadrupoles transverse to the mean flow. The shear- noise term just consists of 
the sum of three independent quadrupoles. 

Equations (47) and (48) show that the self noise depends on the four correlation 
functions R^, R 22 , R i 2 > anc * R 23 - e( l ua ^ on (53) shows that the shear noise de- 
pends on the three correlation functions Rj^, R 2 2 > and ^ 33 - However, it follows from 
equation (40) that 


2 




R 22 R 33 


)df 


(54) 


Hence, the noise is completely determined in terms of four correlations of the turbu- 
lence . 
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ISOTROPIC TURBULENCE 


The expression for the sound intensity can be further simplified if we assume that 
the turbulence is locally isotropic in the moving reference frame. That is, we assume 
that R-j is an isotropic tensor. In this case, the correlation tensor (40) reduces to 

Ky(y,f,T) =A^ j + Bfiy (55) 

where A and B are functions only of £ and r (see ref. 10). 

Upon defining the correlation- length scale L(y , r) by 


L(y , t) = 



df 


(56) 


where P is the rms turbulent velocity at 
appendix C that 


y , we find from the results obtained in 



for i, j = 1, 2, 3 


and 


S 



(57) 


Equation (52) now becomes 


u* y) se 


r 

8ffC 0 ' ’ |_7 


6 -icot g T l dr + G 




(58) 


*It is shown in ref. 4 that the correlation tensor for isotropic turbulence must 
always be an even function of r. Hence, if we began with isotropic turbulence, it would 
not be necessary to make any assumption about the time dependence of R^. 
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In order to simplify the shear-noise term, we use the results of appendix D, which 
show that 


J U'U ,, R n df = J U'U"R 33 df 


and 


J U'U"R 22 df = 0 


If we define the correlation- length scale ./(y,t) by 


Af, t)^ 4/ U’U”R n dr 

Uu 2 


where U denotes the velocity at y , then equation (53) becomes 


(59) 




4P ° u2 “ 2 fck jT e-“"G hJ dr + G 13 £ dr 

c L 


(60) 


Thus, once the Green's function G of the particular geometry has been found, equa- 
tions (58) and (60) give the spectral density of the intensity of the sound field due to a 
unit volume of turbulence located at the point of y as a function of the turbulence 
correlation- length scales L and J!.. 

CALCULATION OF FAR -FIELD GREEN'S FUNCTION 
General Case 

Since we are only interested in the Green's function G(x |y ) when the observation 
point x is in the far field, considerable labor could be saved in calculating the Green’s 
function if the far-field Green's function were calculated directly instead of the general 
Green’s function being calculated and then expanded for x in the far field. This can be 
accomplished by utilizing the reciprocity principle (ref. 13) which states that the Green's 
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function with the observation point at x and the source point at y is equal to the 
Green's function with the observation point at y and the source point at x . That is, 

G(x | y ) = G(y |x ) (61) 

Hence, it follows from equations (11a) and (lib) that 

(V 2 + k 2 )G = -6(x-y) (62) 

n • VG = 0 (on solid boundary) (63) 

where 6(x - y ) = 6(y - x ) and the differentiations are with respect to the variable y 
instead of x . We now put 


G(x |y) = g(x |y) + * (x |y) 

where 


(64) 


g(x|y).__L_e 1 *l?-y| ( 65 ) 

4ff|x-y| 

is the free- space Green's function. Upon substituting this into equations (62) and (63), 
we obtain 


(V 2 + k 2 ) *= 0 

n • = -n • Vg (on solid boundary) 

Upon expanding |x - y | for larger values of |x | we obtain 

|x - y | = r - y + 0(r”*) as r-» 
r 


(66) 

(67) 


where we have put 


r = x 
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It is now convenient to express the location of the source point in terms of the 
spherical coordinates shown in figure 3. Then 

x ' y = yj sin 9 cos ;// + y ^ sin 8 sin i// + yg cos 8 


and 


g(x | y ) ~ — e lkr g (8, \p | y 

4nr 


(68a) 


where 

g Q (0, \}/ |y ) s exp^-ik(yj sin 8 cos if/ + Y 2 si* 1 ^ sin ^ + y$ cos 0)J (68b) 

represents a plane wave propagating in the (-x )-direction. Substituting this into equa- 
tion (67) gives 

n • = - -J— e* kr n • Vg (on solid boundary) (69) 

4 nr 0 

It can now be seen that equations (66) and (69) possess a solution of the form 

*(x|y) =-i_e ikr * (0, ^ |y ) 

4nr 

where 4> 0 is a solution of the boundary- value problem 

(V 2 + k 2 )* o = 0 (70) 

h • V»k = -n • Vg Q (on solid boundary) (71) 

and since 'k represents an outgoing wave at infinity, this is the desired solution to the 
problem. Hence, we can write 

G(x |y) =J_ e ikr G (0,*|y) (72) 

4 nr 0 
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where 


G o (0, \f/ |y) = g 0 (s, |y) + ^ o (0, ^ |y) (73) 

and g Q is the incident plane wave (68) and 'I' Q is the solution to the boundary- value 
problem (70) and (71). It follows from this that G Q is the solution to the problem of a 
plane wave from infinity incident on the solid boundary. By eliminative L he dependence 
on the variable r, the problem has been simplified. 


Two-Dimensional Geometries 

We shall now consider the case where the solid boundaries do not change their shape 
in the y^-direction. Thus, the unit normal vector n lies in the y^yg-plane and 

n • V = n 1 + n« — — (74) 

3y x 3y 2 


The general Green's function of this boundary- value problem is still three dimensional. 
However, the far-field Green's function can be simply expressed in terms of the solu- 
tion to a two-dimensional boundary- value problem. Thus, it follows from equations (68), 
(70) to (72), (73), and (74) that G Q takes the form 

_ -ikyoCosS 

G o (0,tMy) = e G i(^|y 1> y 2 ) (75) 


where 


Gj(^ yj,y 2 ) 


-ik (y^cosi/z+yosim//) 
e + <p 0 (^|y 1 ,y 2 ) 


(76) 


k Q = k sin e (77) 

and the function cp is determined as the solution to the boundary- value problem 
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-ik (y-. cosi//+y 2 sini// ) 

n • V<p Q = -n • Ve (on solid boundary) (79) 

It can be seen from equations (76) to (79) that Gj is the solution to the problem of a 
plane wave from infinity incident on the solid boundary. The wave propagates parallel 
to the y^-plane and has a wave number k Q . 

SUMMARY OF GENERAL EQUATIONS 

Before using these results to work out a particular example, we shall summarize 
the general equations obtained so far. 

The spectral density of the intensity of sound at the far- field point x emitted from 
a unit volume at the point y can be written as 

y* ly) = yxjy>se + y*ly>sh (51) 

self noise shear noise 


Axisymmetric Turbulence 


Self noise: 




e" ia)T G T S dr 


‘ ( G 22 G 33 + G 33 G 22^23 " ^22^ dr 

+ G « £ *' 1 " Tg 3«u dT ] (62) 

Shear noise: 

y ? |y»sh~ foG ij fr^f™**** <53) 

71 
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Isotropic Turbulence 



and the far-field Green’s function G is determined in terms of the incident plane 
problems by 

wave 

General case: 

G(x |y) = J_ e lkr fg (0,i//|y) + (0,^|y)l 

4nr L J 

(80) 


g Q s exp£-ik(y^ sin 9 cos \p + yg sin 9 sin i// + y^ cos 0)J 

(68b) 


(V 2 + k 2 )* 0 = 0 

(70) 


n • = -n ■ Vg Q (on solid boundary) 

(71) 
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Two-dimensional geometries: 


n 


„ 0 3 

V = ni + n 9 

1 3yi 0 y 2 


(74) 


G(x|y) 


t ik(r-y~cos0) -ik (y,cosj//+y 2 sini//) 

— =- e 3 e 0 1 * +^ 0 (^|y 1 .y 2 ) 

4/rr L . 


(81) 


k Q = k sin 8 


3 2 + — + k 2 |<? o = 0 


,ay2 a y 2 


-ik (y 1 cos^/+y 2 sin^/ ) 

n • V(p = -n ■ Ve ° 1 (on solid boundary) 


(77) 

(78) 


(79) 


For both axisymmetric and isotropic turbulence, the shear- noise term is the sum of 
contributions from independent quadrupoles with no coupling between them. The self- 
noise term, however, has this property only in the case of isotropic turbulence since the 
second term in the expression for the intensity for axisymmetric turbulence is a coupling 
term between the two longitudinal quadrupoles transverse to the flow. However, in the 
case where the geometry is two dimensional, it is possible to express this term also as 
a sum of independent quadrupoles. For in this case, equation (81) and the definition of 
G.j show that 

G 33 = -k 2 cos 2 0 G 


and the differential equation (78) shows that 


'll 


°22 * k o G 


0 


Hence, 


G u = tan 2 e G 33 - G 22 (82) 

If this is multiplied by the complex conjugate of the corresponding equation for the time- 
displaced Gjp we find that 
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1 


G 33 G 22 


G 22 G 33 


= tan ^6 


G 33 G 33 


tan^0 


fr r T 
^ 22^22 


' g ii g il 


And this shows that the coupling term in equation (52) can be expressed in terms of inde- 
pendent quadrupoles. 


A PARTICULAR EXAMPLE 

We shall now apply the results obtained in the preceding sections to calculate the 
spectral density of the intensity of the sound emitted from isotropic convected turbulence 
confined between the two semi- infinite parallel plates shown in figure 4. In this case the 
boundary condition (79) becomes 

3 cp -ik y,cos^+ik bsini// 

= ik Q sin ^/e at y 2 = ±b (83) 

3y 2 


The problem of a plane wave incident on two semi- infinite parallel planes subject to the 
boundary condition (83) has been solved by using the Weiner- Hopf technique and is pre- 
sented in reference 14. It is shown there that for |y 2 | < b (see eq. 3.25a in ref. 14) 


<P 


o 



(S_ cosh yy 2 


-yb-iayj 

+ D_ sinh yy 2 )e 1 dor 


(84) 


where the contour integral in the complex a - o + i/3 plane is carried out along any line 
parallel to the a -axis with /3 in the range 

- 9mk Q < /3 < cos \p 9m k Q (85) 

and it is assumed for purposes of obtaining the solution that k Q has been extended into 
the complex plane so that it has a positive imaginary part. At the end of the calculation 
the imaginary part can be set equal to zero. In equation (84), 


y(«) = 




( 86 ) 
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where the branch cuts are chosen as indicated in figure 1. 1 of reference 14. Also, 


D 


k Q sin i {/ cos (k Q b sin 

^2n (k Q + k Q cos ;// ) 1//2 K + (k Q cos i//)(a - k 0 ) 1 ' /2 K_(o')(a - k Q cos \p) 


( 87 ) 


S =■ 


ik Q sin \p sin (k Q b sin ) 


K ± (a) = e 


Vjjff b(k Q + k Q cos i//)L + (k Q cos \p)(a - k 0 )L_(a)(o! - k Q cos i//) 


(88) 


IT fi-^-d/i)] 

n=l ^ 


1/2 ^ iiofb ^/i /ft\ 

Tiffb n-(l/2) e (89) 


L,(a) = e 


*Xo( 


OO 


(90) 


b =2L 

nn 


J n-(l/2) 


( n ' i ) 1 


T t (a)=15yta4±i? 

' V “o / 


where the imaginary part of the logarithm lies between - n and n . 


X X («) = - — (l “ C 0 + In -2- 1 + - ab 


2bk 0 / 2 


X 2 («)=-^(l-C 0 + ln^-J + l«b 
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and C Q = 0. 5772 . . . is Euler's constant. In addition, K ± and L ± satisfy the rela- 
tions 


K (a)K_(a) = e _yb cosh yb 

L + (q;)L_(q!) = (yb)”^e“ yb sinh yb 

at least within the strip (85) of the complex a = a + i/3 plane, 
the appropriate branch cuts shows that 

y(k Q cos \p) = -ik 0 sin ip 

Introducing this relation in equation (91) yields 


(91) 


Using equation (86) with 


(92) 


ik bsim// 

K + (k Q cos ip )K_ (k Q cos \p) = e cos (k Q b sin xp) 




L + (k Q cos \p )L_ (k Q cos xp) = 


ik Q bsini// 
e 

k Q b sin \p 


sin (k Q b sin \p) 


(93) 


Upon substituting equations (88) and (89) into equation (84) and using equations (92) and 
(93) to simplify the results, we obtain 


? o =e 


-ik Q bsin\£/ 



L_(k Q cos t//)(k Q - k Q cos \p) 
L_(or)(Qf - k Q )(Q! - k Q cos \p) 


-yb-iay 1 
e 1 


cosh yyg da 


2/r 


/**>+ u 

•'- 00 +i/3 


K_(k Q cos \p)( k Q 


j yj yj 

K_(a)(a - k 0 )^ // ^(a - k Q cos xp) 




(94) 


Now we are interested in the situation where the source point y is within the chan- 
nel. This corresponds to y^ < 0. In this case, the contour integral can be closed in 
the upper half plane. It can be seen by inspection that the branch cuts in the integrand 
cancel and that the integrand has only poles in the upper half plane. Hence, it can be 
evaluated by the method of residues to obtain 
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= " e 


-ik 0 (yi cos \i/+y2 s i n ^ ) -ikgbsim// 


4 - e 


/ ik o y l\ L - (k o cos ^ 
' / L (k) 


ir>7r i ^n^ b n" y l ^ 

L_(k cos s^)(k - k cos i//)e e y 2 

+ - > COS 


f£ 

n=l 


h -KKK - V^n - k o cos ^ 


n 


oo 

{£ 


K_(k Q cos 4/ )(k Q - k Q cos ip) x e 


l/2 e i [n- ( 1/2 )] (1/2 )[ b n- ( 1/2 f y l] 


sin ■ 


^2 


“ j' K lfn-(l/2)] b n-(l/2)[^n-(l/2) “ k o] 1/2 [ /3 n- (1/2) ’ k o cos *] b n-(l/2) 


where 


r 


2 2 9 kb 

i -* / — k for v > 

b 2 


= < 


kb 


O 2 2 

'k -il- for v ^ — 

O n 


b 


7T 


v = n or v=n — 
> 2 

n = 1,2,3, . . . 


and k Q can now be taken as real. The terms K T J/3^ ( 1 / 2 )] an< * (0 n ) denote 
K_|j3 n _^y2)J an< * L_(/3 n ), respectively, with the factor n = v omitted in the infinite 
product: 


K :[?n-(l/2>] 


1 iQfb n- (1/2) K -» 


L ib n-(l/2) 


a - A 


n_ (l/2)J 


a_/3 n-(l/2) 


L :w n > - 


, iab n L_(ar) 
— e 

ib n " - P n 


a=B 


n 
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Upon using these results in equation (81), we find that the far-field Green’s function 
becomes 


■, ik(r-y,>cos0-bsin0sini^) 

G(x | y ) = — — e J Q{*,e | yi ,y 2 ) 

4nr 


where we have put 


Q(^,o\y 1 ,y 2 ) = e ° yi 




L_(k Q cos 4 /) 
L -< k o> 


(95) 


S 1 J) " i ( " i/3 n )Jj n C0S 7 ' y [- i/3 n-(l/2)] Jl n sin b ~ 

1 \ L-j l ,J V(i/2) 

n=l n=l 


for j = 0, 1, 2 


(96) 


where is the derivative of with respect to y^ and 

(- l) n K (k cos ij/)(k - k cos Yl ] 

j _ " U DO 

K -[^n-(l/2)] b n-(l/2)[^n-(l/2) " k o] 1 ^ 2 [ 3 n-(l/2) " k o cos 

i^n(Vyi) 


V 


(97) 


J n s 


(- 1 ) n L_ (k Q cos 4 /){ k Q - k Q cos ^/)e 

- k o^n - k o cos *> 


These results can now be used to calculate the various quadrupole source terms which 
appear in equations (52), (53), (58), and (59). Thus, 

, ik(r-y^cos0-bsin0sim//) 

G ij(x|y)=— e ^ij(^,0\y v y 2 ) (98) 
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where Q .j - and in view of equation (82) 

Q= Me ° yi + 

Q n = - k o M - e lk ° Yl + s i 2) 

#12 = S 2 1) 

Q = -ik cos #|jik 0 M_e 0 * + 

^22 = tan2 ^ ^?33 ~ Q\\ 

<y?23 = -ik cos 0 S 2°^ 

$33 = -k 2 cos^M.e'* 0 ’' 1 + S<°>] 

and we have put 



J 


(99) 







“ i/3 n-(l/2)] j T _ y 2 

— — I COS 

V(l/2) b n-(l/2) 


for j = 0, 1 

( 100 ) 


We shall now restrict our attention to isotropic turbulence. In order to evaluate the 
integrals in equations (58) and (59) to predict the self and shear noise, it is necessary to 
know the dependence on r of the turbulence correlation lengths L and J. Chu 
(ref. 11) has shown that for the mixing region of a jet 3 L/3t and 3 J'/dr could be 
fit fairly well with functions of the form c^ sech CgT cos c^t (see fig. 37 in ref. 11). 

We shall, for lack of better information about the turbulence, assume that relations of 
this form also hold for the turbulence between the plates. 2 Thus, following Chu we put 


2 Actually, it turns out that the final results are quite insensitive to the details of the 
models used for L and I. 
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( 101 ) 


where 


sech wpV cos u4 2 V 

0T 4 



sech u^r cos a>*^T 


a 


> 



/(IV) 

o 



The time dependence in G and the is associated with the variable y^, which 
only occurs in exponentials. Hence, the integrals over r which must be performed in 
equations (58) and (60) can be found by specializing the constant a in the relation 


r- 


(w-U a)r 


L dr = 


(a; 


, f°° i(w-U a)r 

vi: 


3 4 L 

3t 4 


dr 


L (iv) /*°° i(ci>-U a)r (2) m 

I e sech u>£ 't cos t dT 

« - U P a) 4 J mC0 


t, rL 


(IV) 


2cjJ 2) (o) - U c a) 4 


sech 


3T feu - U a + a)* 2 '] n\oj - U 

_L £ U + sech -1 : 


c a ~ 


(2)1 


2 a) 


( 2 ) 


2a,< 2 > 


and in a similar relation for 2. 

We shall now suppose that the turbulent eddy is far enough from the end of the duct 
so that the contribution to the integrals in equations (58) and (59) are principally from 
the parts of the Green’s function for the interior of the duct. Thus, for the case where 
the eddy is moving toward the end of the duct, we suppose that the point y is far enough 
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from the end so that the eddy at y will decay before it reaches the end of the duct. Us- 
ing values of the eddy lifetime typical of jets, we estimate that this will occur if the 
center of the eddy is about 1/2 to 1/4 of the duct width from the end. When these results 
together with equations (44) and (95) to (100) are used in equations (58) and (59), we find 


V* y> se 


q 4^- t l ™ a (t 


87TCp(477r)^ 


(102) 


l J* >sh = 




nc^(4nr) 2 


where for m = 1,2 and 


£< 2 > = M_D( 2 >e' lk ° yi + < 2 >S<°> 


£<“> = -k 2 M_D (m) e" lk ° yi + (m) S <2) for m = 1,2 

q 8 - (i) 4 1) 

£ j™* = -ik cos ej"-ik 0 M_D (m 'e ° Vl + (ml s ( j 1 for m = 1,: 

,2) ^33 - [2) Qn 


= - ik cos 0 (2) S^ 0) 


(103) 


(104) 


- -k 2 cos 2 e 
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Here we have put 
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and 


M c 



(107) 


is the convection Mach number of the eddy. 

The spectral density of the intensity of the self and shear noise can now be computed 
by substituting equations (89), (90), (95) to (100), and (104) to (106) into equations (102) 
and (103). 

The expressions for I w (x | y ) ge and I (xly)^ become particularly simple when 


37 



\f/ = 0 (i.e. , in the plane of the duct). Thus, it follows from equation (97) that I = J n = 0. 

Equations (95), (96), (100), and (105), therefore, show that all the and are 

equal to zero and M_ = 1. It can now be seen from equations (77), (99) and (102) to (104) that 


and 




ifW 


irc®(4jrr) 2 


at xp = 0 


(108) 
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k 4 D 


^sin 2 * 


at 4> = 0 


(109) 


DISCUSSION 

Emission of Sound From Monopole Source of Unit Strength 


Before considering the sound emission from a region of turbulence within the duct, 
it is useful to consider the sound emission from a monopole source placed in various 
locations within the channel. It is shown in reference 15 (ch. 7) that the intensity of 
sound in the far field from a harmonic monopole source of unit strength can be expressed 
in terms of the far-field Green's function G by 


mon 


AqCqI °| 2 

2tt 


For the two-dimensional channel shown in figure 4, the far-field Green's function is 
given by equations (95), (89), (90), (96), and (97). It can be seen from these equations 
that | G | depends on 6 and k only in the combination 


k sin 6 


Hence, the number of curves necessary for presenting the results can be minimized by 
plotting 
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b 2 sin 2 0 I 
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2 


as a function of i {/ for various values of the wavelength parameter 

A _ 2tt 

b sin 6 kb sin 6 

and various values of the source location y^/b and y 2 /b. The geometrical relation of 
the angles 6 and \p to the duct walls is illustrated in figure 5. Holding 6 and r 
fixed and varying \p corresponds to traversing the circle C in the observation plane p 
parallel to the yjyg- plane. When 8 = w/2, the plane p is the y^y 2 ~ plane. 

The directivity patterns are shown in figures 6 to 8. These figures are polar plots 
2 2 

of (4nrb) sin 8 ^qj/pCq as a function of the angle ^ for various positions of the 
source within the duct. In order to keep the scales of these figures reasonable, the fig- 
ures are divided into three groups. The first group, figures 6(a) to (g), presents the 
curves for wavelength parameters of 8, 6, and 4. The second group, figures 7(a) to (m), 
corresponds to wavelength parameters of 2, 3/2, and 1. And the third set of curves, 
figures 8(a) to (h), are for wavelength parameters of 3/4 and 1/2. The reason that the 
intensities at = n and ^ = -w are not equal is that the duct walls extend to minus 
infinity and thereby separate these points. 

It can be seen from figure 6 that for the long wavelengths the patterns are fairly 
symmetrical, with the sound being radiated fairly uniformly jn all directions. The 
patterns become more uniform as the wavelength increases. However, it can be seen 
from equation (95) that only the first term makes any contribution to the far-field Green’s 
function G for very long wavelengths because the infinite sum becomes very small. 
But this term is the solution for the far-field radiation pattern of a plane wave traveling 
down the duct from yj = (ref. 14). Notice also that this term is independent of the 
vertical location y 2 of the source within the channel, and the intensity I mon is inde- 
pendent of the coordinate yj. Hence, the far-field radiation pattern is independent of the 
source location within the duct for long wavelengths. 

The terms in the infinite sum S; 0 ' represent different modes which can propagate 
within the duct. The first sum in equation (96) represents symmetric modes, and the 
second sum represents antisymmetric modes. As A decreases (i.e., k Q increases), 
the exponentials which enter these terms through equation (97) change from real to com- 
plex (by means of /3 y ) for successively higher order terms in the sums. This corre- 
sponds to successively higher order modes propagating within the duct. Thus, as the 
wavelength becomes shorter, more modes can propagate and, as shown in the figures, 
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the radiation pattern becomes more nonuniform. As each new mode comes into play, 
the pattern changes markedly. 

It can be seen from figure 6, that the first mode in the sum becomes significant 
when X/(b sin 6 ) « 4 if the source is off axis (yg * 0). However, if the source is on 
axis, it does not become significant until X/(b sin 6 ) has decreased approximately to 2. 
This is caused by the fact that only the symmetrical modes can propagate when the 
source is on axis and these modes switch on at shorter wavelengths than the antisym- 
metrical modes. This effect also causes the radiation pattern to be more complex when 
the source is off axis (yg # 0). Thus, as shown in figures 7 and 8, large lobes develop as 
the source is moved from the axis towards the wall. It can be seen that as the wavelength 
becomes shorter, the pattern becomes more complicated. However, the basic outlines 
of these patterns tend to follow the line of sight from the source as the geometric acous- 
tics limit (X « b) is approached. Actually, the wavelengths considered here are too 
long to evidence strong geometric acoustics effects. 


Emission of Sound From a Unit Volume of Turbulence Located 
at a Point Within the Duct and Proposed Experiment 

In the preceding subsection the emission from a monopole source within a duct was 
discussed. In this section we shall discuss the results for the sound emission for con- 
vected locally isotropic turbulence within the rectangular duct. This configuration was 
chosen because of its resemblance to the augmentor wing ejector flap shown in figure 1. 
However, it must be emphasized that the results given in this subsection do not by any 
means represent the total noise emitted from an augmentor wing. They give only the 
sound emitted from a unit volume of turbulence located in the mixing region within the 
augmentor ejector. In order to find the complete noise pattern, it is necessary to inte- 
grate these results over the entire interior of the flap. In addition to this, there are 
other possible noise sources which can contribute to the noise pattern. 

However, it is possible to propose an experiment which corresponds closely to the 
present analysis. This experiment consists of placing a small intense jet at various 
locations within a large wide duct and measuring the narrow-band intensity spectra in 
the far field. Although it is anticipated that the lobes of the observed directivity pat- 
terns will not be nearly as sharp as those predicted by the present analysis (due to the 
fact that the source is not localized and that the microphones do not measure a narrow 
enough frequency range), it is still felt that reasonable agreement should be obtained. 

Discussion of equations. - It has been shown that the spectral density of the inten- 
sity of the sound in the far field emitted from a point y within the duct can be ex- 
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pressed as the sum of two terms. The first of these, I w (x |y) , is called the self 

noise and represents the noise emitted by the action of the turbulent eddies with them- 
selves. This term is given by equation (102). It is proportional to the sum of the in- 
tensities of all six independent quadrupoles corrected for convection effects by the mean 
flow. (Notice that the first term is really a simple source, but it is by virtue of equa- 
tion (99) proportional to the longitudinal quadrupole with axis in the y^-direction. ) As 
in the case of jet noise, this term is also proportional to the fourth- order two- point tur- 
bulence velocity correlations which are accounted for by the factor 

The second term contributing to the far-field intensity is the shear noise I (x |y ) ^ 
which represents the sound due to the action of the turbulence on the mean shear. It is 
given by equation (103). It can be seen that it is proportional to the two longitudinal 
quadrupoles suitably modified by convection effects in the directions perpendicular to 
the velocity gradient (y^ and y^ directions). This term is also proportional to the 
second-order two- point turbulence velocity correlations, as well as the mean-square 
velocity. The turbulence velocity correlations are accounted for by the term 

In the plane 6 - n/2 (shown in fig. 5), a number of the quadrupoles make no con- 
tribution to either the shear or self noise. The reason is the cos 6 dependence of 
these quadrupoles, as shown in equations (99) and (104). Thus, the shear noise is sim- 
ply proportional to the longitudinal y^-direction quadrupole. An inspection of the equa- 
tions for 6 = tt/2 shows that this quadrupole also tends to dominate the self noise at 
low frequencies. 

It is also of interest to have expressions for the intensities in the y^y^-plane, which 
is perpendicular to 6 = 7r/2-plane, as shown in figure 9. This corresponds to setting 
\p equal to 0 and letting 6 vary. The expressions for I w (x|y) ge and I w (x|y) g j i for 
the case where xp = 0 are given by equations (108) and (109), respectively. These 

equations take on particularly simple forms. The directivity behavior of the self noise 

( 2 ) 

is strictly through the term D v , whereas the directivity of the shear noise is through 
the terms D^sin^O. The terms and are given by the first of equations (106). 
They contain the same convective amplification factor (1 - M sin 6)~ 4 that appears in 

c m (2) 

jet noise (ref. 7). In fact, the angular dependence of the terms D v ' and D v is caused 
strictly by the convection of the sound source. The results can also be shown to be con- 
sistent with those obtained by Ffowcs Williams in reference 3. 

Organization of plots and choice of parameters. - The shear- and self-noise terms 


I co (? |y ) se 87 Ic o (47rr)2 
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and 


y*|y > sh i,c o (4,r)2 

2p 0 U 2 ?if> 

are plotted as a function of the azimuthal angle i// for 6 = tt/2 in figures 10 to 16 and 
for 6 = 7i/4 in figure 17. They are plotted as a function of 0 for \p = 0 in figures 18 
to 22. These variations correspond to varying the observation point in the far field in 
the three planes shown in figures 5 and 9. These directivity patterns are influenced by 
the wave numbers of the turbulence which emits the sound. In the present model the 
turbulence spectra are characterized by the parameters k^ m ^and k* m ^ for m = 1 or 
2 (compare eq. (101) where the corresponding frequencies are defined). For the self- 
noise m equals 2, whereas for the shear noise m equals 1. 

In reference 11, it is found that the turbulence data for round jets were best fitted 
by taking 

t^k} 1 ) *M* 
b Q k^ « 2.3 M* 
b^ * 2 M* 
b 0 ki 2) » 5.4 M* 

where M* is the jet- exit Mach number and b Q is the jet radius. For lack of better 
information, these relations will also be used in the present case. We shall suppose 
that M* is the exit Mach number of the flow into the duct. For example, for the aug- 
mentor wing shown in figure 1, M* might correspond to the velocity of the jet entering 
the augmentor flap at point 1. Again, as in a jet, we shall relate the eddy convection 
Mach number to M + by assuming that 
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Figures 1 and 4 show that it is appropriate to take as positive for modeling the 
sound coming from the rear of the augmentor wing exit, whereas it is appropriate to 
take M c as negative for modeling the sound coming from the opening in the front of the 
flap (entrance). Therefore, results for both these cases are presented. For the sound 
from the front of the flap, it would seem appropriate to choose b Q equal to the width of 
the passage at point 1, which is approximately 1/4 b. For the sound coming from the 
rear of the flap, a value closer to b would seem appropriate. Therefore, results are 
presented for both the cases b Q /b = 1/4 and b Q /b = 1. 

Discussion of results . - Results are presented for the 6 = tt/2, 6 - 7r/4, and 
ip = 0 planes. 

6 = tt/ 2 plane: Figures 10 to 16 are plots of the far- field self and shear noise in the 
6 = tt/2 plane. The first two parts of these figures are polar plots of the self noise in 
two different frequency ranges, and the last two are polar plots of the shear noise for 
two different frequency ranges. Comparing the first and last parts of these figures 
shows that the shear noise is always elongated in the direction of the duct axis (y^- 
direction), whereas the self noise tends to radiate uniformly or laterally. This down- 
stream beaming of the shear noise is caused by the fact that, in the 6 = v/2 plane, only 
the longitudinal quadrupole in the y^- direction makes any contributions, whereas the self 
noise is composed of this plus a number of other longitudinal and lateral quadrupoles. 
These figures also show that at any given frequency, the shear-noise curves have much 
fewer lobes than the self- noise curves. This can also be attributed to the longitudinal 
quadrupole nature of the shear noise since the high directivity of this quadrupole in the 
downstream direction allows it to be beamed out of the duct without interacting with the 
walls. There is, therefore, much less interference due to reflections and diffraction by 
the duct than in the more uniformly directed self noise. 

Figures 10 and 11 are identical except that the direction of the convection velocities 
M is reversed. The self- noise curves with M in the positive direction have much 

L v 

weaker lobes than those with M c in the negative direction. The reason for this is that 
the convective effects tend to beam the self noise in the direction of motion. Thus, when 
is in the positive direction, the sound is beamed out of the duct. This beamed pat- 
tern will then interact to a lesser extent with the channel walls than when M £ is nega- 
tive. However, the already highly directional shear-noise patterns show no significant 
additional beaming due to convection, and their directivity patterns are therefore quite 
similar. Figures 12 and 13 show the effect of decreasing the Mach number M* from 
0.8 to 0.3 (with a corresponding decrease in M c ). Comparing figure 12 with figure 10 
shows that increasing the magnitude of the convection velocity in the negative direction 
causes a marked decrease in intensity with no significant change in the shape of the 
directivity patterns. Comparing figure 13 with figure 11 shows that increasing the con- 
vection velocity in the positive direction causes a slight decrease in intensity with no 
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significant change in the shape of the directivity patterns. 

Most of the curves are drawn for the case where b Q /b = 1/4. Figure 14 shows the 
effect of putting b Q /b = 1. A comparison of figures 11 and 14 shows that the curves for 
b Q /b = 1 exhibit less interference effects than those for b Q /b = 1/ 4. 

Figure 15 shows the effect of moving the turbulent source region away from the cen- 
ter of the duct toward the upper wall (y 2 > 0). Comparing this figure with figure 10 
shows that this has little more effect on the shear noise than beaming it in a downward 
direction. However, the self noise exhibits a complex change in the interference effects 
as the source is moved toward the wall. This is probably related to the fact that the 
antisymmetric mode can propagate within the duct when the source is off center. 

Figure 16 shows the effect of moving the turbulent source point farther into the duct. 
Comparing figures 10 and 16 shows that (for both self and shear noise) this causes very 
little change in either the directivity pattern or the magnitude of the intensity. 

6 = ir/4 plane: Figure 17 is a plot of the far- field intensities for the 9 = tt/ 4 plane. 
This plane is (as shown in fig. 5) parallel to the 6 = ir/2 plane but shifted over from the 
centerline of the duct. Upon comparing figures 10 and 17, we see that the shear-noise 
curves on these two planes are similar. However, the figures show that the self-noise 
intensity patterns exhibit weaker interference effects in the 6 = 7i/4 plane than in the 
6 = tt/ 2 plane. 

if/ = 0 plane: Figures 18 to 22 are plots of the far-field, self and shear noise in the 
ip = 0 plane, which corresponds to the y^,- plane shown in figure 9. These plots are 
arranged in the same way as those for the 6 = it/2 plane. Thus, the first two parts of 
these figures are polar plots of the self noise in two different frequency ranges, and the 
last two are polar plots of the shear noise for two different frequency ranges. They also 
correspond to the same values of M c and M* as those for the 6 = u/2 planes. These 
figures are essentially plots of equations (108) and (109). It can be seen from these equa- 
tions that the far-field intensities are independent of the source location within the duct. 
Hence, unlike the previous cases, these coordinates no longer appear as parameters in 
the plots. In addition, the quadrupoles in this plane take on a simple form and, hence, 
cannot exhibit the interference effects which occurred in all other planes. Thus, the 
curves in figures 18 to 22 are much smoother than those in the previous figures. 

These figures also show that the self noise is more intense in the y^direction (i.e. , 
the direction along the duct axis) than in the direction perpendicular to this when the 
convective Mach number M is positive, and that it is less intense along the y^direction 
than in the perpendicular direction when M c is negative. The first two parts of fig- 
ures 19 and 21 show that this elongation is much more pronounced at high convective 
Mach numbers than at low Mach numbers. This is essentially the downstream beaming 
effect caused by convective amplification, which occurs in jets. In fact, the order of 
magnitude of the effect is about the same in the present case as it is for jets (ref. 8). 

On the other hand, figures 10 to 16 show that the shapes of the directivity patterns in the 
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6 = 7t/2 plane are much more dominated by diffraction effects than by convection. Thus, 
the far-field intensity pattern behaves in much the same way as for jet noise in the plane 
of the duct and behaves in a much different fashion in the plane perpendicular to the duct. 

The variation of the intensity with frequency in the if/ = 0 plane is much the same 
as in the 6 = n/2 plane. This is not surprising since the points at 8 - w/2 in figures 18 
and 21 correspond to the points with if/ = 0 (from the same values of M c , M + , and fre- 
quency) in the figures for the 6 = ti/2 plane. As in the case of jet noise, the results for 
0 = 0 are independent of M c (ref. 8). 


SUMMARY OF RESULTS 

A theoretical description has been given of the sound emitted from a nearly parallel 
turbulent shear flow confined to a region near solid boundaries. The analysis leads to 
formulas which can be used to calculate the directivity of the spectra of the sound 
emitted from a unit volume of turbulence located at an arbitrary point within the flow. 

In order to use these results to calculate the sound emitted from any finite region of 
turbulence, it is necessary to integrate the formulas over the region of turbulence. 

This requires a knowledge of certain turbulence- correlation- length scales at each point 
within the flow. 

The results are applied to obtain the sound emitted from a unit volume of turbulence 
located at an arbitrary point in a turbulent flow confined to the interior of an open-ended 
semi-infinite duct. Numerical results for the spectral density of the far-field intensity 
of the self and shear noise are obtained for this case. It is shown that the directivity 
patterns in the plane of the duct are much the same as for jet noise; whereas those in 
the plane perpendicular to this are considerably different, being dominated more by dif- 
fraction effects than by convection effects. This is especially true at the higher fre- 
quencies. The shear noise in the plane perpendicular to the duct is dominated by the 
longitudinal quadrupole with its axis along the duct. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 11, 1972, 

501-04. 
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APPENDIX A 


REDUCTION OF SELF-NOISE CORRELATION INTEGRALS 
FOR AXISYMMETRIC TURBULENCE 


It follows from equation (41) that 


/ i dr = f[ A \Wl + AB(6 ik^j^ + 8 JZ*i«k> + AC(VliMk + 6 li 6 lk^j^> 


+ + 8 u«j«k h + °ik€i€j€z) + D < 8 u 8 ij*k*z 

+ 8 lk 8 lzMj + 8 li 8 lZ*k*j + 6 lk 6 lj^i^ ) + B2s ik 6 j 2 + BC(d ik 6 lj 5 lZ 
* 6 ji 5 li 5 lk* * BD<5 ik 6 lj^ * 5 ik S U^j + 6 ji 6 li 5 k + 6 jt 5 lkh* 


* c S li 6 lj S lk 5 lZ + CD < S H 5 lkV; * 6 li 8 n 5 lk«j + 5 li 6 lj 6 U«k 


Vik 8 u«i>] * 


(Al) 


It can be seen by inspection of this equation that a suitable choice for the coefficient a ik 
in equation (42) is 

a ik - / a2 «?« 2 <f + / (2AD| i + D2 >( 6 ii«k * 6 ik«. 2 ) d f <A2) 


It follows from this equation that 

a ik = a ki 
a 13 = a 12 

But upon using equations (A3) and (A4) we get 


(A3) 

(A4) 
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)G 33 G 33 + a 12 G ii G jj 


+ G 33 G u) + a 23( G 22 G 33 + G 33 G 22) 

+ G 33 G U + G 22 G 33 + G 33 G 22) + (a 23 " a 12 } ( G 22 G 33 
G 33 G 33) + (a 23 ‘ a 12 ) ( G 22 G 33 + G 3 3 G J 2 ) 

+ ^ a 23 ' a 12 ) ( G 22 G 33 * G 33 G 22) 


+ G 33 G 22) 

(A5) 


By using equations (43) and (A5) and the remarks following equation (44), it can now 
be seen by inspection of equation (Al) that 


G ij G kZ M = a 12 G ii G jj + (a 23 “ a 12 } ( G 22 G 33 + G 33 G 22) 

+ G n G'[ 1 ja 11 - a 12 J (2AB{f + 2AC«J + 2D 2 «f + B 2 + 2BC + 4BDi; j + C 2 + 4CD^) dfj 

+ G 12 G 1 2 f (4A 2 £ 2 £ 2 + 2AB^| + 2AB£ 2 + 2AC^ 2 + 2D 2 «| + 8AD^ 2 + 2B 2 + 2BC + 4BD£ j) df 
+ GjgG^y (4A 2 £ 2 £ 2 + 2AB? 2 + 2AB£ 2 + 2AC? 2 + 2D 2 ? 2 + SADt^j + 2B 2 + 2BC + 4BD* j) df 
+ G 22 G2 2 |a22 ' a 12 + f (2AB£ 2 + B 2 ) dfj + G 33 G 33 |a 33 - a 12 + j" (- 2AB^ 3 + B ) d^ j 
+ G 23 Gj 3 / (4A 2 | 2 ^ 3 + 2AB| 2 + 2AB4 2 + 2B 2 ) df 


By carrying out the integrations in polar coordinates, it can be verified that the coeffi- 

T T T 

cients of G^G^ and G^G^ are equal and that the coefficients of ^22^22 and 

T ^ O /*2 7 T 2 

G 33 G 23 are equal (see eq. (A2) and recall that / cos d 1 // = / sin i// dip). 

*4 *7) 


Hence, it follows from equation (A2) that 
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f K ik 8 jZ df = a 12 G ii G jj + (a 23 - a 12>( G 22 G 33 + G 33 G 22) 

[/ 


g u g Ti [J^t + + + 2 AB ^i + 2 AC *i + b2 + 2BC + 4BD ^ + C ' + 4CD?1> ^ ' 3l2 ] 

+ (g u gJ 2 * G 13 c[ 3 )/ (4A h\tl - 2AB t\ + 2ABk\ ♦ 2AC? 2 + 2D 2 * 2 + BAD^l ♦ 2B 2 4 2BC ♦ 4BD£ ,) df 
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(< 


T T 

G 22 G 22 + G 33 G 33/ 


(A6) 


4 G 23 G 2 3 / (4A ^2^3 + 2AB «2 + 2AB ^3 + 2B > d « 


Now it follows from equation (40) that 

Rjl = A 2 ^ + 2AB£j 2 + 2AC£ 2 + 4AD£ 2 + B 2 + 2BC + 4BD£j + C 2 + 4CD£ j + 4D 2 ^ 2 


~2 . 2 t 2 t 2 0Ar , £ t 2 ^ ^ 2.2 

Rl2 = A 5^2 + 2AD^^2 + ^ ^2 

2(R 2 2 + H 11 R 22 ) = 4A 2 | 2 ? 2 + 8AD4^ 2 + 2D 2 4 2 + 2AB£ 2 


+ 2ACij 2 + 2AB4 2 + 2B 2 + 2BC + 4BD£ l 


R 2 2 = A 2 s| + 2AB£ 2 + B 2 



A 2 4 2 4 2 


Hence, it follows from these relations and equation (A2) that equation (A6) can be written 
as 
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G ij G L / *ik^ d ^ 


G ii G jj^*^12 + ( G 22 G 33 + G 33 G 22) f (^23 " ^12) d ^ 

+ G U°Tl/fr ~ ^12) d ^ + 2 ( G 12 G 12 + G 13 G u) ffcl2 
+ ( G 22 G 22 + G 33 G 33)y*(^22 ' ^12) d ^ + 2G 23 G 23 f (*22 


+ *11*22) 

- R 23 ) d£ 


df 


The particular form of the last integral has been obtained by introducing polar coordi- 
nates and carrying out the integration over the angle variable. 
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APPENDIX B 


REDUCTION OF SHEAR -NOISE CORRELATION INTEGRALS 
FOR AXISYMMETRIC TURBULENCE 


It follows from interchanging names of dummy indices in equation (49) that 
GijG \ t f U'lTRj; df = GjjG 1 t f VVHft df * G U G^ f U'U"B df 

♦ (GjjG^ + G U G U'U"D^j df + G n Gn/u'U"C df 

(Bl) 

Upon noting that A is an even function of for i = 1, 2, 3; that U’ and U" are 
functions only of and that any inte e ral of an odd function of ^i for i = 2 ’ 3 1S 

zero, we find that 


Oy«»L «T - GyGT/u'U-A^ df 

And since D is an even function of 4 2 anc * £3 an( * an function of 

(GjjGJi * G U G-)/ U'U"D?j df = 2G u G'[ i yU'U"D| 1 df 

Hence, equation (Bl) becomes 

G l jG T/ U'U”R^ df = G llG T / 1 U’U'HA^i + B + C + 2D^) df 

+ G 12G ^ 2 f U’U'XA^ + B) df + GjgG^ J U'U”(A$3 + B) df 
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We can therefore write from equation (40) that 


°ii G L / or - u'u"Rjj d r 
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APPENDIX C 


REDUCTION OF SELF-NOISE CORRELATION INTEGRALS 
FOR ISOTROPIC TURBULENCE 


It is shown in reference 4 that the requirements of continuity dictate that 

4A + £ — + 11® = 0 

dl; i S4 

where £ = | f | and that this in turn implies that (see ref. 10) there exists a function F 
such that 


A = - 


1 _ 3F 

2£ 


(Cl) 


1 ?F 
B = F + - % — 

2 


Substituting these results into equation (55) shows that 

/«*,. dr f [( f + i ■ {F, )*ik - i x F '] [( F + t 1 A f : 


d£ 


= 5.. 6. . 

lk 


/ / H 4F ') F fl 




I A 2 — 


- 


d? 


where we have put 


F' = 
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Upon introducing spherical coordinates and carrying out the integrations over the angle 
variables, this becomes 


/ ^ = 4 » 6 ik 6 jj f « 2 ( F + i« F ’)( F + | « F ’) d5 


- ~ <Vk* + Vj; + Vk]> 1 “ l 2 (tF') 2 H 


But 


/v (- H H ?r ) « = 4 2 i s FF ' ^ (t F ' )2 ] 


d£ 


F 2 £ 3 




.jrfrFr-i^Y«YJ 0 ?2<SF,)2d{ 


(C2) 


Now we must require that the integral /*R.. df remain finite. Hence, it follows 

J J /•OO 

from equations (55) and (Cl) that this implies that J F£ 2 d£ remain finite. Hence 

l) y 

F must approach zero faster than f 3 as £ - °° . This shows that the integrated term 
in equation (C2) vanishes and we obtain 

/ R ik K j; <f =|[ 5 ik 5 j; + - < 6 ij 5 ki + Vji + £ S 2 (?f') 2 d| 

and since, in particular, 


/r 2 u d f = — /°° ^ 2 (^F*) 2 d£ 

' 11 15 J 0 


it follows that 


/Xk g K <T [66 ik 6 jz ♦ 6 ij5w ♦ d f 
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It now follows from equations (47) and (48) that 


and 



for i, j =1,2,3 
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APPENDIX D 


REDUCTION OF SHEAR -NOISE CORRELATION INTEGRALS 
FOR ISOTROPIC TURBULENCE 


Substituting equation (56) into equation (55) shows that 

/ U’U"R df = / U’U" [f + — £F' - — — F'| df 

J \ 2 2 « / 


(no sum on j) 


Upon introducing cylindrical coordinates with the axial direction along the 42 -axis and 
l? = 4 j + £ 3 , we find that 

J U’U"R 11 df = J U'U”R 33 df 


and 
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where the dependence of F on t has been suppressed. Hence, 

J U , U"R 22 dr = 0 
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APPENDIX E 


SYMBOLS 


coefficient in tensor 
defined by eq. (106) 
constant 

coefficient in eq. (42) 
coefficient in tensor 

one- half of distance separating semi- infinite parallel plates 
coefficient first appearing in eq. (42) 
jet radius 

b/nn, n = 1, 2, 3, . . . 

- - - > n = 1, 2, 3, . . . 

( n ' i> 

coefficient in tensor R^ 

Euler's constant, 0.5772 . . . 
coefficient in eq. (42) 

speed of sound in medium of density p Q and pressure p Q 
constants 

coefficient in tensor R^ 
defined by eq. (106) 
defined by eq. (87) 


F defined by eq. (Cl) 

G(x |y ) Green's function at x due to point source at y 
G jj(x|y) 3 2 G(x|y)/Sy. 0y. 

C*ij(y) Gjj(x|y), where x dependence has been suppressed 
G J(y) G*(y + U c t) 

G^(x j y ) spectral density of Green’s function 
G o (0,^|y) defined by eqs. (72) and (73) 
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Gj(^ ly^J^) defined by eqs. (75) and (76) 

Q defined by eq. (95) 

Cj L defined by eq. (98) 

defined by eqs. (103) and (104), m = 1, 2 
g(x |y) free- space Green's function 

g Q (0,^[y) wave function for plane wave (eq. (68)) 
h typical length associated with distances between boundaries 

I (x) spectral density of far- field intensity at x 

I w (x (y ) spectral density of far- field intensity at x due to unit volume of source 

located at y 

I mon far- field intensity from a monopole source of unit strength 

I n defined by eq. (97) 

i unit vector in y^- direction 

defined by eq. (97) 

K ± (ar) defined by eq. (89) 

K l[^v-(l/2)] K -[^v-(1/2)] with * actor corresponding to n = v omitted 


k 

wave number 


k< m) 

O 

o 

' — - <-*-< 
3 


k o 

k sin 6 


k (m) 

4 m) /e 0 


L 

turbulence correlation length defined by eq. 

(5 0 

L ± («) 

defined by eq. (90) 


l:o v ) 

L (£ v ) with factor corresponding to n = v omitted 

L dv) 

0 

fourth time derivative of L evaluated at r 

= 0 

1 

turbulence correlation length defined by eq. 

(60) 

4 ” 

fourth time derivative of 2 evaluated at r 

- 0 


l typical correlation length associated with turbulent eddies 


M c convection Mach number of eddy 

M Mach number characteristic of turbulence 

e 
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M* 

M _(*) 

n i 

n 

P 

P 

P ij 

p o 

% 

R 

r W*’ T ’ t) ’ 

R ijkz ( y^’ t) 
Ry(y»r, r) 

R i, jk 

R ijkz^^ T) 


jet- exit Mach number 
defined by eq. (95) 

component of unit normal in i th direction 
unit normal 

Fourier transform of P - P Q 

pressure 

pressure tensor 

average pressure 

defined by eq. (47) 

radius of curvature in direction of flow 

fourth- order moving-axis velocity correlation functions defined in 
eq. (28) 

second-order moving-axis velocity correlation function (eq. (34)) 
third- order moving- axis velocity correlation function defined in eq. 
fourth-order moving-axis velocity correlation function defined by eq. 


n(°) 0 ( 1) 1 


ijkZ’ MjkZ 3 

(2) Q + 

ijkZ ,K ijkZ ; 


£ c , 

K ijkZ 


s 

T 


fourth-order velocity correlation functions (eqs. (17) to (20) and 
eq. (24)) 

distance of observation point from origin 


s 

defined by eq. (48) 

s 

defined by eq. (88) 

sg) 

defined by eq. (96) 

c(i) 

S 2 

defined by eq. (100) 

(m)g(j) 

defined by eq. (105) 

(2) s (j) 

b 2 

defined by eq. (105) 


surface 

a large interval of time 


T-. 

ij 


Fourier transform of t- 



T ± (a) 

t 

u(y 2 ) 


it 


U" 


u. 


U, 


u i 


U. 

1 


u 


v,v 2 


V 


v i 


v i 


n 


x l’ x 2’ x 3 

x 


yi,y 2 .y 3 

y 

y' 


y" 

Of 

r 


function of a appearing in expressions for K^a) and L ± (a) 
time 

mean velocity in y^- direction 



convection velocity of eddy 

i-U 

component of turbulent velocity in i n direction 

u i(y’>t) 

UjCy’St + t) 

mean- square turbulent velocity 

volumes for source integration 

magnitude of velocity 

component of velocity in 1 direction 

v i (y’,t) 

v^y’St + t) 

normal velocity component 
components of x 
coordinate of observation point 
components of y 
coordinate of source point 

y +- (ini - v) 

2 1 

-H- 1 ^ 

y + — (irji + T] ) 

2 1 

complex wave number variable 
imaginary part of a 

defined by equation immediately preceding eq. (95) 
defined by eq. (86) 

Kronecker delta 
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S(x - y) 


V 

V 

V 

e 

x 


(i) 


P 

r 

p 

Po 

a 


T 


< ^o ( ^l y l’ y 2 ) 

Xj(a) 

X2(«) 

^ (x | y ) 
* o (0,i//| y) 


U) 




Dirac delta function with argument (x - y) 

magnitude of if 

t th 

component of if in i in direction 

y" - f 
y' - y” 

polar angle measured from y^-direction 
wavelength 

radial coordinate in cylindrical coordinates 

V- U C T 

density 

ambient density 
real part of a 
time translation 
stress tensor 

solution to boundary- value problem given by eqs. (78) and (79) 
function appearing in expression for K ± (a) 
function appearing in expression for L ± (a) 
component of Green's function (eq. (64)) 

factor occurring in (x | y ) which is solution to eqs. (70) and (71) 
azimuthal angle in y^yg- plane measured from y^-axis 
angular frequency 

characteristic angular frequencies appearing in time-dependent part of 
velocity correlation function (eq. (101)) 


Subscripts: 

se self- noise component 

sh shear-noise component 

Superscripts: 


time average 
complex conjugate 
components of vector 


1,2,3 
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Figure 6. - Polar plot of I47rrb sin 9) 2 I m0n /pc 0 as function of <|/ - low frequency range. 
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Figure 6. - Concluded. 
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Figure 7. - Polar plot of (47rrb sin 8) 2 I mon /pc 0 ) as a function of - intermediate-frequency range. 
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Figure 7. - Continued. 
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(b) yj/b = -0.25; y 2 /& s °‘ 75 ’ (dlyj/b « -4.0; y 2 /b s 0.5. 

Figure 8 . - Polar plot of ( 47 itb sin eft^/pcQ as a function of ip - high-frequency range. 
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Figure 8. - Concluded. 
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(c) Shear noise, low-frequency range. (01 shear noise, high-frequency range. 

^M Ur ^n^o " Pattern for spectral density of far-field intensity as function of azimuthal angle q; in 0 » 7t/ 2 plane for convection Mach number of eddy M r *- 0 . 4 ; jet-exit Mach number 

M* = 0 . 8; y^b = - 1 ; y 2 /b- 0 ; and radius ratio b 0 /b- 1/4. c 
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(d> Shear noise, high-frequency range. 
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I^sh ^ 4 ^ 2 l (jJ tttysU 71c o { * m) 

(c) Shear noise, low-frequency range, <d) Shear noise, high-frequency range. 

Figure 13. - Polar plot of directivity pattern for spectral density of far-field intensity as a function of azimuthal angle ip in 9 - v/l plane for convection Mach number of eddy 
jet-exit Mach number M* - 0. 3; yj/b ■ -1; y 2 /b - 0; and radius ratio b 0 /b - 1/4. 
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(a) Self noise, low-frequency range. 
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(b) Self noise, high-frequency range. 

Figure 15. - Polar plot of directivity pattern for spectral density of far-field intensity as function of azimuthal angle in 8 • i\fl plane for convection Mach 
number of eddy M c * -0.4; jet-exit Mach number M* - 0.8; y^/b » -1; y 2 /b - 3/4; and b 0 /b = 1/4. 
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Figure 15. - Concluded. 
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(c) Shear noise, low-frequency range. M Shear noise, high-frequency range. 

Figure 16. - Polar plot of directivity pattern for spectral density ol far-field intensity as a function of azimuthal angle « in 9 ■ ir/2 plane for convection Mach number in eddy M c ■ -0.4; jet-exit Mach number 

M* * 0.8; y j/b « -8; y^/b * 0; and radius ratio * 1/4. 
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Figure 18. - Polar plot of directivity pattern for spect 
number of eddy Mj- s *0. A; jet - exit Mach number 


ral density of far-field intensity as function of polar angle 0 in W=0 plane for convection Mach 
■ 0. 8; and radius ratio bg/b ■ 1/A. 
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Duct 



(a) Self noise, low-frequency range. 



Figure 20. - Polar plot of directivity pattern for spectral density of far-field intensity as function of polar angle 9 in V ■ 0 
plane for convection Mach number of eddy M c * -0. 15; jet-exit Mach number M* 1 0.3; and radius ratio b 0 /b * 1/4, 
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(a) Seif noise, low-frequency range. 


(bl Self noise, high-frequency range. 


Duct 



(cl Shear noise, tow-frequency range. 


Figure 21. - Polar plot of directivity pattern for spectral density of far-field intensity as function of polar angle 0 in 
M c ■ 0. 15; jet-exit Mach number M* - 0. 3; and radius ratio b 0 /b - 1/4. 
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(c) Shear noise, low-frequency range. Shear noise, high-frequency range. 

Figure 22. - Polar plot of directivity pattern for spectral density of far-field intensity as function of polar angle 0 in = 0 plane for convection Mach number in eddy 
M =0.4; jet-exit Mach number M* = 0. 8; and radius ratio b Q /b* 1. 



